(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     10600,        288]
NotebookOptionsPosition[      9856,        260]
NotebookOutlinePosition[     10308,        278]
CellTagsIndexPosition[     10265,        275]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"ThetaD", " ", "=", " ", "46.4551"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"PhiD", " ", "=", " ", "119.4077"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"GammaD", " ", "=", " ", "261.8"}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"GMST", " ", "=", " ", "0.0"}], ";"}], 
  " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Alpha", " ", "=", " ", "240"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Delta", " ", "=", " ", "60.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"LMST", " ", "=", " ", 
   RowBox[{"GMST", " ", "-", " ", "PhiD"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"h", " ", "=", " ", 
   RowBox[{"(", 
    RowBox[{"Alpha", " ", "-", " ", "LMST"}], ")"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Psi", " ", "=", " ", "0"}], ";"}]}], "Input",
 CellChangeTimes->{
  3.4602201864392776`*^9, {3.4602202203765087`*^9, 3.460220228915267*^9}, {
   3.460220303310273*^9, 3.460220344477591*^9}, {3.4602203985177927`*^9, 
   3.4602206849488*^9}, {3.460220811993926*^9, 3.4602210536262417`*^9}, {
   3.460221526202045*^9, 3.4602215415381413`*^9}, {3.460221578492743*^9, 
   3.460221715084841*^9}, {3.460222022045568*^9, 3.460222027580758*^9}, {
   3.460222134382347*^9, 3.460222134866921*^9}, {3.46022382902677*^9, 
   3.4602239092018013`*^9}, {3.460225267912022*^9, 3.460225285782387*^9}, 
   3.460229729349709*^9, {3.460307304601954*^9, 3.460307313049519*^9}, {
   3.4603074282085752`*^9, 3.460307474254854*^9}, {3.4603076303494873`*^9, 
   3.460307649804813*^9}, {3.4603083062926893`*^9, 3.460308326961337*^9}, {
   3.460308435122406*^9, 3.460308474340724*^9}, {3.460308513094954*^9, 
   3.460308518259997*^9}, {3.460404828408684*^9, 3.4604048317202253`*^9}, {
   3.460404994178746*^9, 3.460405002135816*^9}, {3.4604050636791973`*^9, 
   3.46040507825845*^9}, {3.460405722761196*^9, 3.460405724877388*^9}, {
   3.46040582386073*^9, 3.460405857493029*^9}, {3.460405889065001*^9, 
   3.460405936613216*^9}, {3.460405982652748*^9, 3.460405983179892*^9}, {
   3.4604078855223207`*^9, 3.460407888266736*^9}, {3.4604079773862677`*^9, 
   3.460407977751575*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"DegToRad", " ", "=", " ", 
   RowBox[{"Pi", "/", "180"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"ThetaD", " ", "=", " ", 
  RowBox[{"ThetaD", "*", "DegToRad"}]}], "\[IndentingNewLine]", 
 RowBox[{"PhiD", " ", "=", " ", 
  RowBox[{
   RowBox[{"-", "PhiD"}], "*", "DegToRad"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"GammaD", " ", "=", " ", 
   RowBox[{"GammaD", "*", "DegToRad"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Alpha", " ", "=", " ", 
   RowBox[{"Alpha", "*", "DegToRad"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Delta", " ", "=", " ", 
   RowBox[{"Delta", "*", "DegToRad"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"h", " ", "=", " ", 
   RowBox[{"h", "*", "DegToRad"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Psi", " ", "=", " ", 
   RowBox[{"Psi", "*", "DegToRad"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.460307299556074*^9, 3.460307300843107*^9}, {
   3.460307335029076*^9, 3.4603073930348587`*^9}, {3.460307436741077*^9, 
   3.460307457148737*^9}, {3.46030762710186*^9, 3.4603076361474867`*^9}, 
   3.460407780229738*^9, {3.460407882877388*^9, 3.460407920004713*^9}}],

Cell[BoxData["0.8107944493432178`"], "Output",
 CellChangeTimes->{{3.4604078693176622`*^9, 3.460407944625699*^9}, 
   3.460407981106038*^9}],

Cell[BoxData[
 RowBox[{"-", "2.0840575172336333`"}]], "Output",
 CellChangeTimes->{{3.4604078693176622`*^9, 3.460407944625699*^9}, 
   3.460407981114208*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"ArcTan", "[", "0.014151032697933528", "]"}], "*", 
  RowBox[{"180", "/", "Pi"}]}]], "Input",
 CellChangeTimes->{{3.460407795294345*^9, 3.4604078240065536`*^9}}],

Cell[BoxData["0.8107403348488978`"], "Output",
 CellChangeTimes->{{3.460407802684908*^9, 3.460407824562829*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"A", "=", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"-", 
       FractionBox["1", "16"]}], 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "GammaD"}], "]"}], 
      RowBox[{"(", 
       RowBox[{"3", " ", "-", " ", 
        RowBox[{"Cos", "[", 
         RowBox[{"2", " ", "ThetaD"}], "]"}]}], ")"}], 
      RowBox[{"(", 
       RowBox[{"3", " ", "-", " ", 
        RowBox[{"Cos", "[", 
         RowBox[{"2", " ", "Delta"}], "]"}]}], ")"}], 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "h"}], "]"}]}], " ", "-", " ", 
     RowBox[{
      FractionBox["1", "4"], 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "GammaD"}], "]"}], 
      RowBox[{"Sin", "[", "ThetaD", "]"}], 
      RowBox[{"(", 
       RowBox[{"3", " ", "-", " ", 
        RowBox[{"Cos", "[", 
         RowBox[{"2", " ", "Delta"}], "]"}]}], ")"}], 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "h"}], "]"}]}], " ", "-", " ", 
     RowBox[{
      FractionBox["1", "4"], 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "GammaD"}], "]"}], 
      RowBox[{"Sin", "[", 
       RowBox[{"2", "ThetaD"}], "]"}], 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "Delta"}], "]"}], 
      RowBox[{"Cos", "[", "h", "]"}]}], " ", "-", " ", 
     RowBox[{
      FractionBox["1", "2"], 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "GammaD"}], "]"}], 
      RowBox[{"Cos", "[", "ThetaD", "]"}], 
      RowBox[{"Sin", "[", "h", "]"}]}], " ", "-", " ", 
     RowBox[{
      FractionBox["3", "4"], 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "GammaD"}], "]"}], 
      SuperscriptBox[
       RowBox[{"(", 
        RowBox[{"Cos", "[", "ThetaD", "]"}], ")"}], "2"], 
      SuperscriptBox[
       RowBox[{"(", 
        RowBox[{"Cos", "[", "Delta", "]"}], ")"}], "2"]}]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"B", "=", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"-", 
       RowBox[{"Cos", "[", 
        RowBox[{"2", " ", "GammaD"}], "]"}]}], 
      RowBox[{"Sin", "[", "ThetaD", "]"}], " ", 
      RowBox[{"Sin", "[", "Delta", "]"}], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "h"}], "]"}]}], " ", "+", " ", 
     RowBox[{
      FractionBox["1", "4"], 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "GammaD"}], "]"}], " ", 
      RowBox[{"(", 
       RowBox[{"3", " ", "-", " ", 
        RowBox[{"2", " ", 
         RowBox[{"Cos", "[", 
          RowBox[{"2", " ", "ThetaD"}], "]"}]}]}], ")"}], " ", 
      RowBox[{"Sin", "[", "Delta", "]"}], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "h"}], "]"}]}], " ", "-", " ", 
     RowBox[{
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "GammaD"}], "]"}], " ", 
      RowBox[{"Cos", "[", "ThetaD", "]"}], " ", 
      RowBox[{"Cos", "[", "Delta", "]"}], " ", 
      RowBox[{"Cos", "[", "h", "]"}]}], " ", "+", " ", 
     RowBox[{
      FractionBox["1", "2"], 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "GammaD"}], "]"}], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "ThetaD"}], "]"}], 
      RowBox[{"Cos", "[", "Delta", "]"}], " ", 
      RowBox[{"Sin", "[", "h", "]"}]}]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\n", 
 RowBox[{"Fp", "=", " ", 
  RowBox[{
   RowBox[{"A", " ", 
    RowBox[{"Cos", "[", 
     RowBox[{"2", " ", "Psi"}], "]"}]}], " ", "+", " ", 
   RowBox[{"B", " ", 
    RowBox[{"Sin", "[", 
     RowBox[{"2", " ", "Psi"}], "]"}]}]}]}], "\[IndentingNewLine]", 
 RowBox[{"Fc", "=", " ", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "A"}], " ", 
    RowBox[{"Sin", "[", 
     RowBox[{"2", " ", "Psi"}], "]"}]}], " ", "+", " ", 
   RowBox[{"B", " ", 
    RowBox[{"Cos", "[", 
     RowBox[{"2", " ", "Psi"}], "]"}]}]}]}]}], "Input",
 CellChangeTimes->{
  3.4602201864392776`*^9, {3.4602202203765087`*^9, 3.460220228915267*^9}, {
   3.460220303310273*^9, 3.460220344477591*^9}, {3.4602203985177927`*^9, 
   3.4602206849488*^9}, {3.460220811993926*^9, 3.4602210309860163`*^9}, {
   3.460221086634513*^9, 3.4602212103041897`*^9}, {3.4602217581934*^9, 
   3.4602218781805887`*^9}, {3.460221921200897*^9, 3.460221963815404*^9}, {
   3.460222061232535*^9, 3.460222063835287*^9}, {3.460224746550274*^9, 
   3.4602247603990707`*^9}, {3.4602248964472723`*^9, 
   3.4602251707548647`*^9}, {3.460229710636957*^9, 3.4602297542151537`*^9}}],

Cell[BoxData[
 RowBox[{"-", "0.18005343506558352`"}]], "Output",
 CellChangeTimes->{{3.460229711026325*^9, 3.460229756498817*^9}, 
   3.460229792309061*^9, 3.460308540378903*^9, 3.4604048475137453`*^9, 
   3.460405006404138*^9, 3.460405069790276*^9, 3.4604057313925943`*^9, 
   3.460405842322694*^9, {3.4604059003479223`*^9, 3.460405916762815*^9}, 
   3.460405988074955*^9, 3.460407947241716*^9, 3.4604079840936127`*^9}],

Cell[BoxData["0.9278375899650511`"], "Output",
 CellChangeTimes->{{3.460229711026325*^9, 3.460229756498817*^9}, 
   3.460229792309061*^9, 3.460308540378903*^9, 3.4604048475137453`*^9, 
   3.460405006404138*^9, 3.460405069790276*^9, 3.4604057313925943`*^9, 
   3.460405842322694*^9, {3.4604059003479223`*^9, 3.460405916762815*^9}, 
   3.460405988074955*^9, 3.460407947241716*^9, 3.460407984138856*^9}]
}, Open  ]]
},
WindowSize->{1179, 789},
WindowMargins->{{115, Automatic}, {Automatic, 22}},
ShowSelection->True,
Magnification->1.25,
FrontEndVersion->"6.0 for Linux x86 (32-bit) (June 19, 2007)",
StyleDefinitions->FrontEnd`FileName[{"Creative"}, "NaturalColor.nb", 
  CharacterEncoding -> "iso8859-1"]
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[568, 21, 2236, 46, 277, "Input"],
Cell[CellGroupData[{
Cell[2829, 71, 1190, 27, 227, "Input"],
Cell[4022, 100, 140, 2, 51, "Output"],
Cell[4165, 104, 157, 3, 51, "Output"]
}, Open  ]],
Cell[4337, 110, 195, 4, 51, "Input"],
Cell[4535, 116, 112, 1, 51, "Output"],
Cell[CellGroupData[{
Cell[4672, 121, 4342, 121, 499, "Input"],
Cell[9017, 244, 420, 6, 51, "Output"],
Cell[9440, 252, 400, 5, 51, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
